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ABSTRACT 

We propose a new interpretation of the quasar luminosity function (LF), derived from physically motivated 
models of quasar lifetimes and light curves. In our picture, quasars evolve rapidly and their lifetime depends 
on both their instantaneous and peak luminosities. We study this model using simulations of galaxy mergers 
that successfully reproduce a wide range of observed quasar phenomena. With lifetimes inferred from the 
simulations, we deconvolve the observed quasar LF from the distribution of peak luminosities, and show that 
they differ qualitatively, unlike for the simple models of quasar lifetimes used previously. We find that the 
bright end of the LF traces the intrinsic peak quasar activity, but that the faint end consists of quasars which are 
either undergoing exponential growth to much larger masses and higher luminosities, or are in sub-Eddington 
quiescent states going into or coming out of a period of peak activity. The "break" in the LF corresponds 
directly to the maximum in the intrinsic distribution of peak luminosities, which falls off at both brighter and 
fainter luminosities. Our interpretation of the quasar LF provides a physical basis for the nature and slope of 
the faint-end distribution, as well as the location of the break luminosity. 

Subject headings: quasars: general — galaxies: nuclei — galaxies: active — galaxies: evolution — cosmology: 
theory 



1. INTRODUCTION 

The luminosity function (LF) of quasars is fundamental to 
cosmology, bu t even after more than 30 years of study (e.g., 
Schmidtlll968t ISchnudt & Greed fl983l lHartwick & Schade 
1990 1 IWarren et all Il994t iBovleetalJ 120001 lUedaetalJ 
120031) . its relation to the intrinsic luminosities of quasars 
remains poorly understood. Previous work modeling the 
quasar LF has relied on restrictive assumptions about life- 
times and light curves of quasars, imagining, for example, 
that quasars either have universal lifetimes or that they evolve 
exponentially, usually on the galaxy dynamical time or the 
e-folding time for Eddington-limited black hole growth t$ = 
Mbu/M = 4 x 10 8 e r yr for accretion with radiative efficiency 
e r = L/Mc 2 - 0. 1 ( SalpeterH%l. Under these circumstances, 
the distribution of quasars with a given mass or peak luminos- 
ity is trivially related to the observed LF (in the absence of se- 
lection effects), and the two have essentially identical shape. 

Recently, simulations of gal axy mergers incorpo rating 
black hole growth and feedback (Springel et al. 2005b) have 
reproduced the Mbh - a rel ation between black ho le mass and 
galaxy velocity dispersion (Di Matteo et al. 2005) and linked 
quasar activity ( Hopkins et al. 2005a b) to galaxy evolution 
dSpringel et alJl2005al) . The simulations predict more com- 
plicated quasar light curves than have been adopted previ- 
ously. The peak, exponential black hole growth phase is 
determined by the gas supply over timescales ~ 10 s yr and 
shuts down when significant gas is expelled b y feedback. Th e 
light curves have been studied by Hopkins et al. (2005a b), 
who showed that the self-termination process gives observ- 
able lifetimes ~ 10 7 yr for bright optical quasars, in good 
agreement with observations, and yields a large population 
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of obscured sources a s a natural stage of quasar evolution. 
Hopki nset alJ l l2005bl) analyzed simulations over a range of 
galaxy masses and found that the quasar light curves and life- 
times are always qualitatively similar, with both the intrinsic 
and observed quasar lifetimes being strongly decreasing func- 
tions of luminosity, with longer lifetimes at all luminosities 
for higher-mass (higher peak luminosity) systems. 

Here, we use our quasar lifetimes and light curves to decon- 
volve the observed quasar LF from the distribution of peak lu- 
minosities, and find that they differ qualitatively, unlike for the 
trivial light curves or even complex, cosmologically evolving 
light curves or distributions of Eddington ratios that have been 
employed previously, demanding a new interpretation of the 
quasar LF. 

2. MODELING QUASAR LIFETIMES & LIGHT CURVES 

Up to now, theoretical studies of the quasar LF 
have generally employed very simple descriptions of the 
quasar light curve, namely some vari ant of a "feast or 
famin e" or "light bulb" model (e . g.. ISm all & Blandford 
[T99l iKauffmann & Haehiieltl l2000l IWvithe & Loebl l2003t 
Haiman. Ouataert. & Bower 2004), in which quasars have 
only two states: "on" or "off." When "on," quasars have a 
constant luminosity and radiate at a characteristic Eddington 
ratio / = (L/LEdd), generally in the range ~ 0.1 — 1. This 
state is assumed to last for a universal lifetime f@./.B, which 
is the same for all quasars at a given redshift (although some 
models allow for redshift evolution of ?q,lb). The value of 
tQ,LB is an input parameter of the models, generally adopted 
from observations or assumed to be related to the Salpeter 
time or the dynamical time of the host galaxy. The time spent 
in a logarithmic luminosity interval is then just dt/dlog(L) = 
tQ.LB ln(10)L5(L— lL,Edd)- Although this approach is analyti- 
cally simple, with the trivial light curve 

f{t) = lL Edd Q{t)Q{t Q . LB -t) 1 (1) 
where 9 is the Heaviside step function, it has no strong 
theoretical or observational motivation. For present pur- 
poses, models where quasars live arbitrarily long with 
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slowly evolvin g mean volume emissivity or mean light 
curve (e.g.LSmall & Blandfordll992tlHa iman & Menou 2000; 
Kauffmann & Haehnelt 2000) are equivalent to the "light 
bulb" scenario, as they still assume that quasars observed at a 
luminosity L radiate at that approximately constant luminos- 
ity over some universal lifetime tQ.LB at a particular redshift. 

We also consider a variant of the "light bulb" scenario, 
which we term the "Eddington limited" model, where a black 
hole accretes at a fixed Eddington ratio / from an initial mass 
Mi to a final mass Mj (or equivalently, a final luminosity 
Lf = ILEddiMf)), and then shuts off. This gives exponential 
mass and luminosity growth, with the light curve 

f(t) = lL Edd {M i )e""\ (2) 

where /(f) = for t > InMf/M/. The time spent in any loga- 
rithmic luminosity bin is constant, 

dt/d\og(L) = t s (\n(10)/l) (3) 

for Li < L < Lf. This is true for any exponential light curve 
/(f) oc e ±x l u , such as that of lHaiman & Loebl dl998t) . with 
only the normalization dt /d\og(L) = f* changed, and thus any 
such model will give identical results with correspondingly 
different normalizations. While such a light curve allows the 
black hole mass and luminosity to change, there is no theoret- 
ical expectation that black hole growth is well-approximated 
by a constant Eddington ratio. 

We compare these scenarios to a physically motivated 
model usi ng quasar lifetimes der ived from simulated galaxy 
mergers (lHor)kinsetal.ll2005alhh . The light curves from the 
simulations are complex, generally having periods of rapid 
accretion after "first passage" of the galaxies, followed by 
an extended quiescent period, then a transition to a peak, 
highly luminous quasar phase, and then a dimming as self- 
regulated mechanisms expel gas from th e remnant center afte r 
the black hole reaches a critical mass (Di Matteo et al. 2005). 
In addition, the accretion rate at any time can be variable 
over small timescales ~ Myr, but even with these complex- 
ities, the statistical nature of the l i ght curv e can be described 
by simple forms. Hop kins et al.l (l2005bl) find that the total 
quasar lifetime fg(Z/ > L) above a given luminosity L is well- 
approximated by a truncated power law, with 

t Q (L> >L) = t 9 (L/l0 9 L Q r, (4) 

where f 9 = t Q (L' > 10 9 L Q ) ~ 10 9 yr, over the range 10 9 L Q < 
L < L pea k for a given quasar. L pea i is, as above, determined 
by the final black hole mass, L pea t rj L^ddiMf)- Over a wide 
range of L peak (from - 10 10 - 10 14 L^). lHopkins et al.l (T2005b) 
also find that a is a function of My (or L pea k), given by 
a = ao + a' \ogL pea k, with a = -0.2 (the approximate slope 
of the Eddington-limited case) as an upper limit. This reflects 
the fact that larger quasars have shallower slopes as they spend 
more time at higher luminosities up to some larger peak lumi- 
nosity. The time spent in any logarithmic luminosity interval 
in this range is then simply 

dt/dlog(L)=\a\t 9 (L/l0 9 L Q ) a . (5) 

In our picture, quasars spend far more time at low luminosi- 
ties than at their peaks. This is a consequence of both pure 
1=1 Eddington-limited growth up to a final black hole mass 
as well as a "brightening" and "dimming" out of and into ex- 
tended quiescent phases of black hole growth around the peak 
quasar phase, during which the Eddington ratio changes with 
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FIG. 1 . — A typical quasar light curve (top) and lifetime as a function of 
luminosity (bottom), for a major merger of Milky Way-like galaxies with a fi- 
nal black hole mass Mf = 3 X 10 8 Mq i Hopkins et al. 2005a). The simulation 
lifetime (solid line) is compared to the expected Eddington-limit model life- 
time (dotted), the power law approximated lifetime (dashed), and the mixed 
power law-Eddington model (dash-dot), from the formulae in §2. 



time in the range I ~ 0.0 1 - 1. We also consider a more com- 
plicated fit described in Hopkins et al. (2005b), where quasar 
lifetimes follow the Eddington-limited model at the bright- 
est luminosities, down to some L», below which they are de- 
scribed by a power law. In all these cases, the best-fit break 
luminosity is a nearly constant fraction of the peak luminos- 
ity L* w Q.\AL pea k. These fits have the advantage of being 
accurate at both low and high luminosities, down to lifetimes 
<Myr. In Figure 1, we show an example of the light curve 
from one of our simulations (top) and the relation between 
the lifetime and luminosity and various fits to it (bottom). 

Our model contrasts sharply with both the "light bulb" 
and "Eddington limited" cases. The "light bulb" scenario 
ignores the evolution of the quasar as it accretes, and ne- 
glects the vast majority of its life at luminosities below its 
final peak. The "Eddington limited" model may be a rea- 
sonable approximation during the bright quasar phase near 
the final peak luminosity, but it underestimates the time 
spent at lumino sities L < 0.1 L nm u b y as much as two orders 
of magnitude ( iHopkins et alJl2005bl) . Furthermore, in our 
model, quasars can evolve significantly on short timescales 
(~ 10 Myr), which means that our interpretation yields a sta- 
tistical mix of contributions from different peak luminosities 
and Eddington ratios at any observed luminosity at a given 
redshift, which contrasts strongly with any model in which 
the quasar light curve or emissivity evolves over cosmic time. 
We note that our arguments do not in any way depend on the 
obscuration model described in Hopkins et al. (2005af3), as 
we consider only low enough redshifts that quasar LFs should 
be reasonably complete and therefore only "intrinsic" quasar 
lifetimes are relevant. Incorporating such effects, however, 
does not qualitatively change our conclusions. Moreover, 
the simulations from which we derive our lifetimes cover the 
complete range of intrinsic and observed luminosities of the 
LFs considered here, and show smoothly-changing properties 
in good agreement with the fits described above for all simu- 
lations in this range. 

3. RESULTS 
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FIG. 2. — Deconvolved intrinsic distributions of peak luminosities 
from fitting to the lBovle et alj i2000) LF (diamonds) at redshift z = 0.5. Up- 
per panel shows the fitted LF, lower panel the ri(L„ ea k) distribution for "light 
bulb" models (dashed), "Eddington limited" models (dotted), and luminosity 
dependent lifetimes fitted to power laws (thin) and fitted to combined power 
laws and Eddington growth (thick). 



Given quasar lifetimes as functions of luminosity and peak 
luminosity, the LF is a convolution of the lifetime with the 
intrinsic distribution of sources with a given L pea k, 



$(L) oc 



dlog(L) 



ri(L peak )dlog(L peak ). 



(6) 



Here, n(L peak ) is the rate at which sources in a logarithmic in- 
terval in Lp ea k are created or "activated" per unit volume at 
some redshift, and Q(L) is the observed number density of 
sources per logarithmic interval in L. This formulation im- 
plicitly accounts for the "duty cycle" (the fraction of active 
quasars at a given time), which is proportional to the lifetime 
at a given luminosity. As we do not attempt to model the 
number of halos containing quasars and wish to consider wide 
classes of lifetime models (which may have different normal- 
izations), we do not concern ourselves with the absolute nor- 
malization of Q(L) or n(L pea k). 

For the "light bulb" and "Eddington limited" scenarios re- 
spectively, this convolution and the corresponding deconvolu- 
tion are trivial, giving 
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for the intrinsic L pea k distributions. ISteed & Weinberg! J2003I) 
use a similar method to de-convolve the distribution of 
sources and consider in detail a very wide range of models for 
the quasar light curve and distribution of Eddington ratios, but 
none of the models they consider include the fundamentally 
important feature of the lifetime increasing with decreasing 
luminosity to the lowest observed quasar luminosities, and 
further they do not allow the Eddington ratio distribution to 
depend on final black hole mass. As a result, every example 
in the range of models they consider, as well as those in pre- 
vious works, results in an n{L pea k) distribution qualitatively 
similar to the observed LF, increasing monotonically to lower 
luminosities (in the absence of an arbitrary truncation). How- 
ever, in the case of our luminosity-dependent lifetimes, we 
must solve numerically for n LDL , fitting to a given LF, and 
find a fundamentally different result. 

Figure [2] shows the be st-fit n(L nen k ) distributions for the 
above cases fitted to the iBovle et all ([2000) optical quasar 



LF at redshift z = 0.5, as well as the best-fit LF. For 
luminosity-dependent lifetim es, we tak e a pow er law de- 
pendence on Lpeak from iHopkins et alJ (2005B, with a = 
-0.95 + 0.32 log(L^/10 12 L Q ) for pure truncated power-law 
fits and a = -0.98 + O.321og(L / , efl(t /lO 12 L ) for fits with a 
power-law mapped onto an Eddington-limited fit at near-peak 
luminosities. For ease of comparison with black hole mass 
and Eddington ratios, we rescale all B-band quantitie s to the 
bolom etric luminosity using the corrections of Marconi et al. 
(2004). In all cases, it is possible to reproduce the double- 
power law LF and "break" luminosity L^k quite accurately, 
and the n{L pea k) distributions are similar at luminosities above 
Lbrk- However, below Lb rk , a proper accounting of the lumi- 
nosity dependence of quasar lifetimes results in a radically 
different n(L pea k) distribution, with the shallow slope of the 
LF a direct reflection of the slope in the lifetime vs. luminos- 
ity relation, dominated by quasars at the peak of the n(L pea k) 
distribution, which then determines the "break" luminosity of 
the observed LF, with L(nix>L.max) ~ L\, r k. Note that the steep- 
ness of the ri(Lp ea k) cutoff is an artifact of extrapolating the 
LF to luminosities well below those observed - any strong 
decline below L^k gives an equivalent fit to the actual data. 

It is possible that objects which have already undergone a 
quasar phase may re-activate, or that black holes may already 
be very large when AGN activity first begins. To determine 
whether or not this changes our results, we consider the most 
ex treme possible case, in which we return to the simulations 
of lHopkins et alJ fe005albl) and re-fit the quasar lifetimes but 
ignore all black hole activity prior to the final merger. This 
neglects all early accretion activity as black holes build up to 
a significant mass, as well as weaker, sub-Eddington accretion 
during the merger as gravitational torques from early merger 
stages funnel gas to the galaxy centers. Therefore, we expect 
this case to set a strong lower limit to the lifetime at low lumi- 
nosities, regardless of the mechanism driving quasar activity. 
However, we still find that these lifetimes are well above the 
Eddington-limited model expectations, albeit with shallower 
best fit power law slopes a = -0.77 + 030log{L pea k/10 n LQ) 
and a slightly lower (constant) normalization. 

It has also been argued from observations of stellar black 
hole binaries that a transition between accretion states oc- 
curs at a critical Eddington ratio m = M/M£dd, from ra- 
diatively ineffic ient accretion flows at l ow accretion rates 
(e.g.. lEsin. McClintock. & Naravanll 19971) to radiatively effi- 
cient accretion through a standard Shak ura~& Sunvaevl i ll 9731) 
disk. Although the critical Eddington ratio for supermassive 
black holes is unc ertain, observations of black hole binaries 
(M accaron as well as theoretical extensions of ac- 

cretion models (e.g., Mey er. Liu. & Mever-H ofmeister l20"o"ol) 
suggest m cr it ~ 0.01. We examine whether this can change 
our conclusions by re-fitting quasar lifetimes from our simu- 
lations. Because we assumed a constant radiative efficiency 
L = e r Mc 2 with e r = 0.1, we account for this effect by multi- 
plying the simulation luminosity at all times by an additional 
"efficiency factor" f e ff which depends on the Eddington ra- 
tio I = L/LEdd (given an a priori constant assumed efficiency, 
I = m always), 



feff : 



1 if/ > 0.01 
100/ if l< 0.01. 



(9) 



Thi s choice for t he efficiency factor follows from ADAF mod- 
els (Naravan & Yi 1995) and ensures that the radiative ef- 
ficiency is continuous at the critical Eddington ratio Z cr „ = 
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FIG . 3 . — As Figurel2l lower panel, but with lifetimes recalculated consider- 
ing only times after the onset of peak quasar activity (thick), or decreasing the 
radiative efficiency below a critical Eddington ratio (thin). The dot-dashed 
line uses lifetimes calculated including both effects. "Light bulb" (dashed) 
and "Eddington limited" (dotted) models are shown for comparison. Trian- 
gles show the distribution of black hole masses expected from the distribution 
of early-type velocity dispersions i Sheth et al. 2003 ) and the Mbh ~ cr relation 
fromfTremaine et al. 1 20021). 



0.01. We expect this to set an extreme limit to the impact 
of this transition, as we do not incorporate this effect dy- 
namically in the simulations, which would slow "blowout" 
after the peak quasar phase and increase lifetimes at low lu- 
minosities. We again find that power-law lifetimes at low 
luminosities are much larger than the expectation from the 
Eddington-limited model, with best-fit slopes a = -0.94 + 
0.331og(W/10 12 L o ). 

Figure [3] shows the best-fit n(L pea k) distributions and LFs 
for these cases, in the manner of Figure[2] We also plot as an 
extreme limit the results obtained applying both the radiative 
efficiency correction above and considering only times dur- 
ing and after the final merger. This case gives the shallowest 
slopes, a = -0. 50+0.5 log(L peak /10 12 L Q ). Although the shal- 
lower slopes broaden the n(L pea k) distribution, our results are 
qualitatively unchanged. We have further considered extreme 
cases such as making our original slopes shallower by 3er, a 
factor ~ 2, and find nearly identical results. 

4. CONCLUSIONS 

We consider realistic quasar light curves and life- 
times derived from hydrodynamical simulations of quasar 
activity, and find that this results in a qualitatively dif- 
ferent form and evolution of the distribution in peak 
quasar luminosities n{L pea k) (equivalently, final black 
hole masses or host system properties ) than that derived 
or implied in previou s work s (e.g., Small & Blandford 
T99l lHaiman & Loebl fl99l IHaiman & Menod | 2000[ 
Kauffmann & Haehnelt T l2000t IWvithe & Loebl \20oi. 
Volonteri et alJ 120031 IHaiman. Ouataert. & Bowed 12004 
We describe the resulting new interpretation of the quasar lu- 
minosity function (LF) in which the steep bright-end consists 
of quasars radiating near their Eddington limits and is directly 
related to the distribution of intrinsic peak luminosities (or 
final black hole masses) as has been assumed previously. 
However, the shallow, faint-end of the LF describes black 
holes either growing efficiently in early stages of activity or in 
extended, quiescent states going into or coming out of a peak 



bright quasar phase, with Eddington ratios generally between 
/ ~ 0.01 and 1. These are not objects which were active at 
some earlier cosmic time and have since "shut down", as 
are invoked in some "light bulb" models, but correspond to 
the evolution of quasars in mergers over short timescales 
and critically related (both physically and statistically) to 
the observed bright, Eddington-limited quasar phase. The 
"break" luminosity in the LF corresponds directly to the peak 
in the distribution of intrinsic quasar properties n(L pea i). As 
such, observations of this break and its evolution present a 
powerful new probe of the quasar population. 

We find that the observed quasar LF can be reproduced to 
high accuracy with a model of quasar lifetimes which depends 
on both the observed luminosity and peak luminosity of a 
quasar, for some distributi on of peak lumino sities n{L pea k). 
We demonstrate this for the lBovle et alJ <2000l) LF at z = 0.5 , 
but find identical qual itative results for the M ivaii et"aT1 d2000l) 
and lUeda et alJ J2003I) soft and hard X-ray LFs, over a wide 
range of redshifts z = 0-3. This is not trivial, as it is in 
the case of simpler models of quasar lifetimes (in which the 
n(L pea k) distribution can be directly recovered from the LF). If 
the faint-end slope of the LF is too shallow, power-law lumi- 
nosity dependent lifetimes will be too steep to reproduce it. If 
the faint-end slope is too steep, our lifetime model can repro- 
duce the LF, but without a characteristic qualitative difference 
in behavior from "light bulb" or Eddington-limited models. 
In fact, it is a strong argument in favor of our picture that ob- 
served faint-end slopes lie just within the expected range, for 
any peaked population n(L pea i). 

Our new interpretation of the quasar LF immediately ex- 
plains several other observations. The distribution of veloc- 
ity dispersions in early-type systems (exp ected to have under- 
gone merger-driven AGN activity) from S heth et all (2003) 
turns over and decreases below a w 160kms _1 , which from 
the Mbh - cr relation implies that the black hole mass dis- 
tribution for remants of major mergers of massive galaxies 
should turn over and decrease below Mbh ~ 1O 8 M (L p ea t ~ 
1O 12 L ). Moreover. IHaiman. Ouataert. & BowerH2004l) find 
that black holes with masses less than ~ 10 7 M Q must be 
rare at high redshift to avoid significantly overpredicting the 
counts of bright radio sources. These observations are exactly 
as predicted by our deconvolution of the LF, where at high 
redshift n{L pea k) is shifted to higher peak luminosity (black 
hole mass) and thus fewer low-mass black holes are pro- 
duced. The n{L pea k) distribution based on traditional models 
of quasar lifetimes, however, would necessarily predict that 
the vast majority of quasars, especially at high redshift where 
the number of quasars is larger, have final black hole masses 
well below this limit, with the number density increasing to- 
wards lower black hole masses. The shift of n(L pea k) with 
time in our model can explain the en tire observed range o f the 
z = black hole mass function (e.g., Marco nTet all2004l) . 

The distribution of observed Eddington ratios is also a natu- 
ral consequence of our model, and does not invoke a particular 
probability distribution across sources to match observations. 
Quasars observed around and below the "break" in the LF 
(which dominate the total population) may be accreting at the 
Eddington rate, as they build up early in their lives, or may 
be in quiescent phases going into or out of their peak quasar 
activity. For these quiescent phases, both comparison of the 
lifetime power laws in §3 and direct calcul ation of the life- 
time as a function of Eddington ratio ( Hopki nTet alJ2005albl) 
show that they are dominated by relatively large Eddington 
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ratios I > 0.1. Below I ~ 0.01, quasars at the peak of the 
n{L pea k) distribution will be below limits measured in most 
surveys. Therefore, we expect a distribution of Eddington ra- 
tios concentrated between I ~ 0.1 — 1, with a small tail down 
to I ~ 0.01 and mean Eddingt on ratios increa sing with lumi- 
nosity, exactly as observed ( Vestergaard 2004]). 

The sharp contrast between the intrinsic distribution of peak 
quasar luminosities and the observed LF provides clean and 
elegant tests of our theory. In our interpretation, the bright 
and faint ends of the LF correspond statistically to similar 
mixes of galaxies, but in various stages of evolution. How- 
ever, in all other competing scenarios, the quasar luminos- 
ity is directly related to the mass of the host galaxy. There- 
fore, any observational probe that differentiates quasars based 
on their host galaxy properties such as, for example, the 
dependence of clustering of quasars on luminosity, can be 
used to discriminate our picture from olde r models. Early 
evidence from both quasar-quasar ICroom et alJ 120051) and 
quasar-galaxy (Adelberger& Steidel 2005) correlations sup- 
port this picture and suggest that quasar hosts have a well- 
defined characteristic mass, exactly as predicted by our model 
with a peaked niLp^) and in stark contrast to previous mod- 
els in which n(L pea k) continues to increase below the break in 
the luminosity function. 

Our modeling of realistic quasar lifetimes and the result- 
ing new liiLpeak) distribution and interpretation of the quasar 
luminosity function have a wide range of implications, . The- 
oretical predictions of the active and relic supermassive black 



hole density and mass functions, the contribution of quasars 
to the X-ray and infrared background, the observed Edding- 
ton ratio distribution as a function of luminosity and redshift, 
the evolution of characteristic black hole masses and quasar 
host galaxy masses with redshift, quasar correlation functions 
and bias, the relation between quasar and starburst or lumi- 
nous infrared galaxy luminosity functions, the role of quasars 
in reionization, the population of high-redshift radio sources, 
means to discriminate between pure density and pure lumi- 
nosity evolution of quasar populations at very high redshift, 
and the distribution and evolution of quasar host properies 
all depend on this quantity and must be revised in this new 
model. Semi-analytical models and cosmological simulations 
have focused on reproducing n(L pea k) distributions qualita- 
tively similar to the observed LF based on idealized models 
of the quasar lifetime, but should instead attempt to repro- 
duce the fundamentally different n(L pea k) distribution implied 
by our realistic, physically motived models of the quasar light 
curve and luminosity-dependent quasar lifetimes. 
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